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Abstract 

Floating  ice  impact  forces  are  of  concern  where  structures  ore  built  in  ice- 
susceptible  waters.  Bridge  piers,  ice  control  structures  or  icebreakers  ore  o  few 
examples  where  the  ability  to  predict  the  expected  ice  impact  force  wouid  be  of 
great  help  in  the  design  process.  Experiments  were  performed  to  determine  the 
response  of  a  floating  ice  beam  to  a  verticaily  applied  force.  The  data  were  used 
to  calibrate  a  numericai  finite  element  model  of  the  floating  ice.  The  ice  was 
characterized  as  a  linear  elastic  material  in  the  numericai  anaiysis,  and  the 
calibration  data  were  used  to  assess  this  assumption  as  well  as  to  develop  o 
fluid  influence  coefficient  matrix  to  simulate  the  dynamic  influence  of  the  fluid 
beneath  the  ice  beam.  Finally,  a  scale  model  study  wos  performed  to  determine 
actual  impact  forces  generotsd  by  a  floating  Ice  beam  ogainst  a  45°  sloped 
structure.  The  numerical  model  developed  was  then  compared  to  the  actual 
data.  The  numerical  model  does  well  at  predicting  impad  forces  for  all  the 
beams  at  low  velocity  and  the  force  from  the  thicker  ice  beams  at  all  velocities. 
Both  the  numerical  and  experimental  forces  show  the  same  trends  and  appear 
to  level  off  and  approach  a  constant  value  with  increasing  beam  length.  The 
discrepancies  between  numerical  predictions  and  experimental  results  are 
thought  to  be  caused  by  damage  in  the  experimental  ice  beams  which  is  not 
qccounted  for  in  the  numerical  model. 


For  conversion  of  SI  metric  units  to  U  .S./Br1tish  customary  units  of  measurement 
consult  ASTM  Standard  E380-89a,  Standard  Pmcticetor  Use  of  the  International 
System  of  Units,  published  by  the  American  Society  for  Testing  and  Materials, 
1916  Race  St.,  Philadelphia,  Pa.  19103. 
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Experimental  and  Numerical  Analysis  of 
Floating  Ice  Beam  Impact  Forces  Against  A  Sloped  Structure 

BARRY  A.  COUTERMARSH 


INTRODUCTION 

The  force  generated  by  a  floating  ice  floe  im¬ 
pacting  against  a  structure  is  of  interest  in  many 
engineering  situations.  Icebreakers,  bridge  piers 
in  cold  regiot\s,  and  ice  control  and  offshore  struc¬ 
tures  are  just  a  few  of  the  structures  that  must 
withstand  ice  impacts.  It  is  difficult,  if  not  impos¬ 
sible,  to  analytically  determine  these  forces  be¬ 
cause  of  the  complex  nature  of  the  ice-structure 
interaction  involving  the  supporting  fluid,  as  well 
as  the  difficulty  in  characterizing  the  material  re¬ 
sponse  of  the  ice. 

The  initial  impetus  for  this  work  was  to  deter¬ 
mine  the  maximum  force  that  could  occur  during 
an  ice  impact  against  the  U.S.  Army  floating  tacti¬ 
cal  "Ribbon  Bridge."  The  bridge  comes  in  sections 
that  are  deployed  individually  and  then  connected 
in  the  water  to  form  a  continuous  bridge.  Each 
section  has  two  rectangular  roadway  pontoons 
and  two  bow  pontoons  with  45°  sloped  faces  as 
shown  in  Figure  1.  During  its  normal  use  in  a  river 
where  floating  ice  might  be  present,  the  sloped 
faces  of  the  bow  pontoons  would  receive  any 
impact. 

This  report  describes  a  model  study  of  ice  im¬ 
pacting  a  45°  sloped  surface.  The  forces  generated 
by  the  impacts  are  analyzed  using  a  general  finite 
element  technique  with  the  ice  characterized  as  a 
linear  elastic  material.  The  fluid/ice  interaction  is 
handled  by  a  linear  fluid  inertial  coupling. 

The  study  was  performed  in  essentially  three 
parts.  Experiments  were  performed  to  measure 
the  response  of  a  floating  ice  beam  to  a  verticaUy 
applied  force  at  one  end  of  the  beam.  The  data 
from  these  experiments  were  used  to  calibrate  a 
numerical  finite  element  model  of  the  floating  ice. 
It  was  used  to  help  assess  the  assumption  of  a 
linear  elastic  material  designation  for  the  ice  and  to 
develop  a  fluid  influence  coefficient  matrix  that 
simulated  the  dynamic  influence  of  the  fluid  be¬ 
neath  the  floating  ice  beam.  This  step  is  described 
in  the  beam  response  section  below.  In  the  next 
step,  a  model  study  was  performed  to  gather  data 
on  the  actual  forces  generated  by  a  floating  ice 
beam.  The  impact  data  were  then  analyzed  using 


a  finite  element  numerical  model.  These  last  two 
steps  form  the  main  body  of  the  work  and  are 
described  in  detail  below. 


BACKGROUND 

The  force  generated  by  an  ice  floe  impacting 
against  the  ribbon  bridge  is  difficult  to  determine 
analytically  due  to  the  unconstrained  motion  of 
both  the  floe  and  bridge. 

The  ribbon  bridge  is  a  series  of  articulated  sec¬ 
tions  consisting  of  four  floating  pontoons  attached 
by  hinges  and  latches  (Fig.  1).  The  bow  pontoons 
offer  a  45°  sloped  surface  for  an  ice  floe  to  impact 
against.  This  surface  is  free  to  respond  to  the  floe's 
kinetic  energy  in  a  number  of  ways: 

1.  The  surface  itself  can  be  indented  or  tom. 

2.  The  bow  pontoon  can  displace  due  to  the  free 
movement  in  the  hinges  and  latches  holding  the 
pontoons  together  as  well  as  through  the  slack  in 
the  cormectors  between  bridge  sections. 

3.  Entire  bridge  sections  can  translate  due  to  the 
slack  in  the  anchoring  system  securing  the  bridge 
in  the  river. 

4.  The  pontoon  can  displace  due  to  the  hydro- 
dynamic  response  of  the  floating  bridge  to  balance 
the  impact  force. 

In  addition  to  the  bridge's  response,  the  impact 
energy  can  be  dissipated  in  the  ice  through: 

1.  Crushing  and  fracturing  of  the  ice  during 
impact. 

2.  The  irreversible  hydrodynamic  response  of 
the  fluid  foimdation. 

3.  Horizontal  rotation  of  the  floe. 

4.  The  internal  strain  energy  respoi\se  of  the  ice. 

Only  the  portion  of  the  impact  energy  that  is 

actually  transferred  to  the  bridge  is  capable  of 
damaging  the  bridge,  either  through  indenting  or 
puncturing  of  the  pontoon  or  failure  of  the  cormec¬ 
tors.  The  numerous  ways  that  the  floe's  energy  can 
be  dissipated  make  it  difficult  to  analytically  quan¬ 
tify  the  amount  that  would  be  trairsmitted  to  the 
bridge  versus  that  which  goes  into  the  other  mecha¬ 
nisms  mentioned  above.  It  should  also  be  pointed 
out  that  since  material  failure  is  of  concern,  it  is 
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Figure  1.  Profile  ofU.  S.  Army  fioating  ribbon  bridge. 


necessary  to  determine  the  maximum  impact  force 
as  opposed  to  the  impulse  force.  This  distinction 
means  that  the  total  amount  of  force  overall  is  not 
important  but  rather  the  maximum  peak  that  oc¬ 
curs  during  the  impact. 


THEORY 

Initially  the  problem  was  looked  at  from  a  rigid 
body  d)mamics  point  of  view  with  the  forces  de¬ 
scribed  by  Newton's  laws.  From  Newton's  first 
law,  the  ice  is  ideally  considered  to  be  in  equilib¬ 
rium  as  it  floats  on  ^e  water  and  moves  with  the 
current  at  a  constant  translational  velocity  and  no 
rotational  motion.  Its  velocity  is  taken  to  be  the 
same  as  the  water  and  it  therefore  has  no  relative 
motion  to  it.  Using  a  reference  axis  on  shore,  I 
considered  the  bridge  to  be  a  rigid  body  at  rest  in 
the  middle  of  the  waterway  with  no  motion. 

Newton's  third  law  states  that  for  any  force 
interaction  between  two  bodies  the  second  body 
always  exerts  an  equal  and  opposite  force  to  the 
first.  The  force  that  is  measu^  at  the  model  is 
therefore  identical  to  the  force  experienced  by  the 
ice  at  the  instant  of  impact.  This  force  is  the  same 
regardless  of  how  the  impact  is  created;  i.e.,  it  is 
irrelevant  if  the  ice  moved  into  the  bridge  or  the 
bridge  moved  into  the  ice. 

The  force  created  during  this  impact  is  deter¬ 
mined  by  the  mass  of  the  bodies  involved,  their 
initial  velocities,  the  acceleration  that  occurs  dur¬ 
ing  the  impact  and  their  velocities  after  the  impact. 

After  the  impact,  if  the  solution  of  the  lineiu' 
velocity  vector  v  and  angular  velocity  vector  oo  of 
the  two  impacting  bodies  were  possible,  then  the 
forces  associated  with  the  impulse  of  the  impact 
could  be  determined  by  summing 

I  F  dt  =  {miVj2  + mb^,i}-{miVil-mbVbl) 

(1) 


for  the  total  linear  momentum  and 


_  _  _  _ 

I  T  dt  =  (ij  (0i2  +  lb  t0b2)  -  (fi  ®il  -  fb  ®bl) 

Jtl  (2) 

for  the  total  angular  momentum.  Since  the  ice  and 
bridge  are  free  to  respond  in  three  dimensions  the 
determination  of  their  velocities  after  the  impact 
would  be  quite  difficult. 

The  above  also  ignores  the  contribution  of  the 
water.  The  fluid  in  contact  with  the  rear  portion  of 
the  floe  will  decelerate  and  displace  as  itstrikes  the 
decelerating  block.  The  fluid  beneath  the  block  is 
also  accelerated  as  the  block  accelerates  vertically 
into  the  fluid.  These  processes  use  energy  and 
would  need  to  be  determined  for  an  accurate 
energy  accoimting. 

Furthermore,  if  frictional  heating,  fracturing  or 
plastic  deformation  occurs  to  either  body,  some 
portion  of  the  energy  involved  in  the  collision 
would  be  used  in  these  processes  and  would  be 
difficult  to  find  analytically. 

Ausefulrelationfordescribingimpulsiveevents 
is  the  coefficient  of  restitution  e,  whi^  depends  on 
the  shape  and  material  properties  of  the  impacting 
bodies.  The  coefficient  of  restitution  is  de^ed  as 
die  ratio  of  the  impulse  force  during  restoration  to 
the  impulse  during  deformation  for  each  body,  or 


where  F,  is  the  contact  force  during  restoration,  F^ 
is  the  contact  force  during  deformation,  f  is  the 
total  time  of  contact  between  the  bodies,  and  td  is 
die  total  time  of  deformation. 

While  the  coefficient  of  restitution  does  not 
detail  the  local  or  internal  behavior  of  the  colliding 
bodies,  it  may  be  able  to  represent  the  bdiavior 
when  determined  either  aiialytically  or  experi¬ 
mentally.  It  is,  however,  dependent  upon  velocity, 
as  well  as  the  material  and  body  configurations. 
For  perfecdy  elastic  bodies  with  no  energy  loss,  e = 
1,  and  for  perfecdy  inelastic  bodies  or  plastic  im¬ 
pact,  e  =  0  and  the  energy  loss  is  at  a  maximum. 

All  of  the  above  problems  make  an  experimen- 
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tal  determination  of  impact  force  desirable.  They 
also  point  out  the  desirability  of  an  accurate  nu¬ 
merical  technique  to  model  impacts. 

MODEL  SIMILITUDE 

The  modeling  procedure  assumed  that  gravity 
forces  were  of  primary  importance.  Thus  to  main¬ 
tain  dynamic  similarity,  the  same  ratio  of  inertial 
force  to  gravitational  force  from  the  model  to  the 
prototype  was  necessary.  If  the  Froude  number, 
defined  as: 


where  V  is  velocity,  g  is  the  gravitational  constant 
and  L  is  some  convenient  length,  is  a  constant 
between  the  model  and  prototype,  then  dynamic 
similarity  is  achieved  (Sharp  1981). 

The  pertinent  similarity  equations  for  the  study 
were  then 

V’m=V'p{LVLp)l/2  (5) 

and  for  the  forces 

fm=  Fp(lrn/Lpf  (6) 

where  F  is  the  force.  The  subscripts  m  and  p  denote 
model  and  prototype  (full-scale)  respectively. 

With  a  1  /5th  scale  model  study  the  above  equa¬ 
tions  reduce  to 

V„=  0.4472  Vp  (7) 

and 

Fn,  =  0.008  Fp.  (8) 

Since  the  prototype  forces  needed  to  be  deter¬ 
mined,  the  above  equation  was  rearranged  as 

Fp=125Fn,.  (9) 

EXPERIMENTAL  PROCEDURE 

The  experiments  were  performed  in  CRREL's 
Ice  Engineering  Facility  test  basin.  The  basin  is  36 
m  long  X  9  m  wide  x  2.6  m  deep  and  the  room  is 
capable  of  being  cooled  down  to  -20°C. 

In  the  study  some  of  the  variables  present  in  full 
.scale  were  removed  to  cut  down  on  the  degrees  of 


freedom  for  the  problem.  The  bridge  was  assumed 
to  be  a  rigid  body,  with  no  internal  flexure  and  no 
external  freedom  to  move  vertically,  or  to  rotate 
about  a  horizontal  axis  perpendicular  to  the  ice 
motion.  The  ice  floes  were  configxured  to  impact 
squarely  against  the  model,  thus  minimizing  en- 
ergy  going  into  horizontal  rotation  of  the  floe. 
These  precautions  assured  that  the  measured  im¬ 
pact  force  would  be  conservative;  i.e.,  it  would 
model  the  highest  force  that  is  expected  to  occur 
naturally. 

A  heavily  reinforced,  1.37-m-wide  wooden 
model  of  a  bow  pontoon  and  a  portion  of  a  road¬ 
way  pontoon  were  constructed.  A  thin  steel  skin 
was  added  to  protect  the  face  of  the  pontoon 
against  ice  abrasion.  The  bow  pontoon  provided 
the  impact  surface  and  the  adjacent  roadway  por¬ 
tion  of  the  pontoon  helped  to  direct  the  fluid  flow 
under  the  bow,  thus  avoiding  any  turbulence  that 
might  not  be  present  in  full  scale  (Fig.  2a).  This 
configuration  allowed  us  to  determine  impact 
forces  against  the  pontoon  but  would  not  model 
forces  at  the  connectors  between  bridge  sections  or 
pontoons. 

The  top  of  the  model  was  rigidly  moimted,  via 
a  force  transducer,  imder  a  moving  carriage  over 
the  test  basin.  The  transducer  measured  forces  in 
three  directions  by  sampling  at  a  constant  rate  of 
500  samples/second.  The  transducer  analog  out¬ 
put  was  converted  to  a  digital  signal  and  stored  on 
a  floppy  disk.  This  allowed  us  to  move  the  bridge 
into  the  ice  and  create  an  impulsive  force  equal  to 
the  naturally  occurring  situation  where  the  ice  is 
propelled  into  the  bridge.  Figure  2a  shows  a  pro¬ 
file  view  of  the  test  setup  with  the  X,Y  and  Z 
directions  indicated. 

The  floes  were  cut  from  freshwater  ice  in  all 
combinations  of  four  lengths  and  thicknesses,  with 
a  constant  width  of  1 .37  m,  as  shown  in  Appendix 
B.  The  sizes  were  chosen  to  span  a  range  that  could 
be  expected  in  field  conditions.  The  model  study 
represented  prototype  floes  of  6.85-m  width,  with 
lengths  varying  from  1.52  to  15.24  m  and  thick¬ 
nesses  from  130  to  760  mm. 

Only  one  size  floe  width  was  used  because  the 
force  was  assumed  to  be  linearly  proportional  to 
width,  as  discussed  later.  The  purpose  of  varying 
the  lengths  and  thicknesses  was  to  determine  any 
effect,  such  as  that  of  ice  flexure  or  rotational  under 
turning,  other  than  the  simple  increase  in  mass  of 
the  larger  floes.  The  impact  tests  were  performed 
at  four  discrete  model  velocities  in  the  expected 
range  of  prototype  water  velocities  encountered  in 
the  field.  The  combinations  of  the  above  variables 


3 


a.  Profile  of  bridge  model,  transducer  and  ice  floe. 


Basin  Walls 


b.  Top  view  of  basin  showing  the  basin  carriage,  ribbon  bridge  model  and  ice  floe  in  the  clear-cut 
slot  in  the  ice. 


Figure  2.  Test  setup.  The  model  was  driven  at  a  constant  velocity  into  the  floating  ice  beam.  In 
the  free  body  diagram,  forces  are  shown  at  the  interface  and  at  the  transducer. 
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resulted  in  64  tests,  in  one  complete  set.  Four  repli¬ 
cates  of  each  test  were  planned,  which  resulted  in 
a  total  of 256  tests.  Some  replicates  and  one  test  were 
not  performed,  so  the  actiial  number  was  246. 

TTie  model  ice  sheets  were  normally  grown  by 
cooling  the  room  air  down  below  freezing  until  the 
desired  t>'.ickness  of  ice  was  obtained.  The  excep¬ 
tion  to  this  procedure  was  the  152-mm-thick  ice, 
which  was  obtained  by  stacking  and  freezing  two 
76-mm-thick  ice  sheets  together. 

A  slot  slightly  wider  than  the  model  bridge  was 
cut  in  the  ice  sheet,  providing  a  clear  path  for  the 
model  to  move  in.  A  floe  used  in  a  particular  test 
was  placed  in  the  slot  and  aligned  parallel  to  the 
adjacent  solid  ice  surface  but  still  free  to  displace 
vertically  and  horizontally  in  the  direction  of  im¬ 
pact.  Two  light  A^eight  sticks  were  placed  on  both 
the  floe  and  the  ice  sheet  to  hold  the  orientation  of 
the  floe  before  impact,  as  shown  in  Figure  2b. 


PRELIMINARY  TESTS 

Preliminary  tests  were  performed  to  determine 
the  repeatability  of  the  experimental  results,  to  test 
the  assumption  that  impact  force  varies  linearly 
with  floe  width,  and  to  test  the  effect  of  bridge 
displacement  upon  impact  force. 

Tests  run  to  determine  the  repeatability  of  the 
experiments  indicated  a  wide  scatter  in  the  results. 
Some  of  the  floes  were  not  impacting  the  bridge 
squarely  and  as  a  consequence  some  amount  of 
energy  was  lost  to  floe  rotation  during  these  eccen¬ 
tric  impacts.  The  repeatability  of  the  tests  was 
improved  by  designing  the  floes  with  a  narrow 
nose.  This  resulted  in  a  head-on  impact  and  gave 
reasonably  consistent  impact  loads. 

The  details  of  the  preliminary  work  are  given  in 
Appendix  A. 

Although  some  discrepancy  was  noted  in  the 
data,  an  assumption  of  linearly  varying  force  with 
floe  width  appears  to  be  valid  from  the  results  of 
preliminary  tests. 

Bridge  displacementhad  an  effect  upon  impact 
force,  with  the  largest  variation  being  at  the  maxi- 
miun  displacement  tested.  At  this  displacement, 
the  average  force  in  the  longitudinal  direction  had 
a  9%  drop,  while  the  average  Z-force  showed  a 
31%  average  increase  over  the  lower  two  displace¬ 
ments  tested.  The  majority  of  tests  were  run  at  this 
maximum  displacement  of  the  model  bridge,  to 
have  a  conservative  value  for  the  vertical  forces 
while  the  longitudinal  force  values  are  perhaps 
10%  low  in  the  worst  case. 


NUMERICAL  PROCEDURE 

A  finite  element  analysis  was  performed  using 
a  combination  of  ABAQUS  software  (Hibbitt, 
Carlson  and  Sorenson  Inc.  1989)  for  the  solid  me¬ 
chanics  (ice  and  bridge)  portion  and  a  fluid  influ¬ 
ence  coefficient  matrix  to  account  for  the  influence 
of  the  fluid  inertia  (McGilvary  et  al.  1990).  The  ice 
beam  was  modeled  using  beam  elements  which 
use  Euler-Bemoulli  beam  theory;  i.e.,  any  traris- 
verse  shear  deformation  is  ignored,  with  a  cubic 
interpolation  of  displacement.  Thenumerical  beam 
was  split  up  into  30  elements  with  a  variable  nodal 
spacing  along  its  length.  A  dense  nodal  spacing 
was  used  in  the  impacted  end  of  the  beeun,  gradu¬ 
ally  becoming  less  dense  by  a  factor  of  5%  per  node 
towards  the  opposite  end  of  the  be<un. 

The  bridge  was  modeled  as  a  massless, 
45°-sloped,  rigid  surface  using  an  element  that 
allows  vertical  and  horizontal  displacement  but 
no  rotation.  Variable  displacement  stiffness  is  in¬ 
corporated  by  attaching  massless  springs  to  the 
surface,  as  shown  schematically  in  Figure  3.  The 
spring  stiffness  was  calculated  by  utilizing  the 
measured  free  vibration  frequency  of  the  bridge 
model,  an  estimate  of  the  effective  mass  of  the 
bridge,  and  the  natural  frequency  equation: 

f=y|:/2n  (10) 

where  /  =  the  oscillation  frequency 
I:  =  the  spring  stiffness 

m  =  the  mass  of  the  bridge. 

An  average  oscillation  frequency  of  13.5  Hz  was 
measured  and,  with  an  estimated  mass  of  34  kg  for 
the  bridge  model,  gave  a  value  of 250,000  N/m  for 
the  stiffness  of  the  surface  support  in  the  majority 
of  the  numerical  runs.  The  ice  elastic  modulus  was 
2.0  GPa  and  it  was  assumed  to  be  isotropic.  A 
series  of  numerical  runs  using  152-mm-thick  ice 
were  performed  to  test  the  effect  of  structure  stiff¬ 
ness  and  ice  elastic  modulus  on  impact  force.  In 
these,  the  spring  stifhiess  was  doubled  to  500,000 
N/m,  to  study  structure  stiBhess,  and  the  ice 
elastic  modulus  was  increased  to  10  GPa,  which 
corresponds  to  a  hard  freshwater  ice  (Ashton  1986) . 

The  static  portion  of  the  fluid  pressure  beneath 
the  ice  is  detoed  as  a  distributed  elastic  foimda- 
tion  (Wiidder  foimdation)  by  a  series  of  springs  at 
each  beam  node.  These  springs  have  a  nonlinear 
spring  stiffriess,  shown  in  Figure  4,  to  simulate  the 
changing  static  fluid  pressure  on  the  ice  (buoy¬ 
ancy  force)  when  it  displaces  from  equilibrium. 

The  dynamic  influence  of  the  fluid  is  calculated 
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Figure  3.  The  finite  element  model  of  the  bridge  and  ice  beam.  U  is 
horizontal  displacement  and  W  is  vertical  displacement.  The  fluid 
foundation  is  represented  by  massless  springs. 


in  a  separate  fiiute  element  procedtire.  In  this,  the 
fluid  domain  of  the  flume  is  discretized  into  finite 
elements  with  an  area  corresponding  to  the  ice 
beam  on  top.  Through  the  wave  equation 

dt^ 


where  c  is  the  speed  of  sound  and  is  the  acoustic 
pressure.  By  ignoring  the  first  term  as  very  small, 
the  acoustic  presstu%  that  is  caused  by  a  unit 
acceleration  of  each  succeeding  ice  node,  in  turn, 
while  the  others  are  held  at  zero,  is  calculated  for 
each  node  in  the  fluid  grid.  The  acoustic  pressure 
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W,  Ice  Deflection 
a.  Static  foundation  reaction  force. 
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W,  Ice  Deflection 

b.  Foundation  modulus  mth  respect  to  displacement. 

Figure  4.  Plots  of  static foundation  reaction  force  and foundation  modulus  with 
respect  to  displacement,  for  an  ice  node,  pj  is  ice  density,  Py,  is  water  density,  h 
is  ice  thidcness  and  g  is  the  graxntational  constant. 
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values  at  the  fluid/ice  interface  are  then  multi¬ 
plied  by  the  area  of  the  ice  over  which  they  act, 
giving  an  inertial  force  at  the  interface.  This  iner¬ 
tial  force  is  then  divided  by  the  unit  acceleration  to 
give  a  fluid  added  mass,  at  each  ice/water  inter¬ 
face  node.  The  resulting  value  of  added  mass  for 
each  node  is  the  fluid  influence  coefficient  matrix. 
This  fluid  influence  coefficient  matrix  is  then  added 
to  that  of  the  ice  beam  in  the  numerical  model, 
where  its  multiplication  with  the  ice  beam  accel¬ 
erations  yields  the  fluid  inertial  forces  acting  on 
the  ice.  More  detail  on  the  fluid  inertial  formula¬ 
tion  can  be  found  in  (McGilvary  1989). 


BEAM  RESPONSE 

To  assess  the  ability  of  the  model  to  simulate  an 
impulsively  loaded,  floating  beam,  a  series  of  tests 
were  done  to  measure  the  response  of  a  floating  ice 
beam  to  a  vertical  impulsive  force  at  one  end.  The 
beams  were  either  54  or  60  mm  thick  and  305  mm 
wide.  The  lengths  varied  from  2  up  to  30  m  and 
typically  had  two  free  ends.  The  characteristic 
length  of  the  ice  is  a  way  of  expressing  the  magni¬ 
tude  of  ice  stiffness  to  foimdation  stiffiiess.  It  was 
determined  by  first  measuring  the  displacement 
caused  by  placement  of  a  1-kg  dead  weight  on  a 
representative  ice  beam  and  using  the  following 
formula  (Hetenyi  1946): 


^~2KB 

where  w  is  the  end  deflection  of  the  beam  caused 
by  the  load  P  at  the  end,  K  is  the  foimdation 
modulus,  in  this  case  the  specific  weight  of  water 
(Ashton  1986),  and  B  is  beam  width. 

The  elastic  modulus,  E,  was  calculated  from  this 
by  applying  the  theory  of  an  elastic  beam  on  an 
elastic  foimdation  with  the  following  formula 
(Hetenyi  1946): 

Lc  =  (£AV12JC)^/^  (13) 

where  E  is  the  effective  elastic  modulus  and  ft  is  the 
ice  thickness. 

The  value  of  was  found  to  be  2.68  m  and  E  to 
bel.72GPa. 

The  floating  ice  beams  were  instrumented  along 
their  length  with  five  linear  variable-displacement 
transducers  (LVDTs).  Typically,  one  LVDT  was 
placed  at  the  impacted  end,  one  about  1  m  back 
from  the  impacted  end,  one  in  the  middle  at  the 
beam,  one  halfway  between  the  middle  and  the  far 
end  and  one  at  the  far  end.  A  vertical  impulsive 
load  was  applied  to  one  end  by  a  wooden  dowel, 
instrumented  with  a  load  cell.  Figure  5  shows 
schematically  a  typical  test  arrangement.  After  the 
short  duration  load  was  applied,  the  free  vibra¬ 
tions  of  the  beam  were  recorded  by  the  LVDTs. 

Figure  6a  shows  a  typical  force-history  trace 


Figure  5.  Schematic  diagram  of  a  floating  ice  beam  instrumented  with  LVDTs.  An 
impulsixx  load  is  applied  at  the  end  by  the  bar  and  measured  by  the  load  cell. 
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b.  Displacement  histories  at  different  points  on  the  beam. 

Figure  6.  Typical  results  from  an  8-m-long  beam  response  experiment. 


created  by  the  applied  vertical  impulsive  load,  and 
Figure  6b  shows  the  measured  response  of  the 
beam  at  several  distances.  The  force  typically  rose 
monotonically  for  0.03  to  0.07  s  to  its  peak  value 
and  then  decreased  monotonically  at  a  somewhat 
slower  rate  to  zero.  Overall  the  force  histories  were 
relatively  symmetrical  about  the  peak  value.  It  can 
be  seen  in  Figure  6b  that  the  impact  end  of  the 


beam  exhibited  the  largest  displacements  over 
time,  as  expected.  The  displacements  further  along 
the  beam  show  a  slight  time  lag  from  the  maxi- 
miun  displacement  at  the  impacted  end  as  well  as 
smaller  displacements  as  the  impulsive  energy  is 
lost  to  the  beam  and  supporting  fluid. 

Figure  7  shows  a  comparison  between  the  ex¬ 
perimental  and  numerical  displacement  data  for  a 
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b.  1.05  m  from  the  impact. 


c.2mfromlhempaa.  i.  4  mfr«m  the  impact. 

Figure  7.  Comparison  of  numerical  and  experimental  beam  displacement  histories  fdra4-m-long  beam. 


4-m-long  beam.  Comparisons  of  theoretical  and 
experimental  responses  of  8-m-long  and  15.25- 
m-long  beams  are  shown  in  Figures  8  and  9,  re¬ 
spectively.  The  model  increasingly  overpredicts 
the  measured  deflections  at  points  farther  from  the 
loaded  end  of  the  beam.  If  one  point  on  the  beam 
is  observed  and  time  allowed  to  progress,  the 
numerical  model  shows  the  best  agreement  in  the 
first  0.5  second.  In  other  words,  the  greatest  dis¬ 
crepancy  between  the  numerical  and  experimen¬ 
tal  data  occurs  at  the  end  of  the  beam  farthest  from 
the  impact  or  after  the  first  0.5  s.  This  is  under¬ 
standable,  because  any  vibrating  floating  body 
will  transfer  kinetic  energy  to  the  fluid  mass  by 
generating  surface  waves,  and  the  numerical  model 
does  not  account  for  this  energy  transfer. 

The  numerical  model  also  assumes  elastic  be¬ 
havior  for  the  ice  and  complete  energy  conserva¬ 
tion.  This  means  that  there  is  no  loss  of  energy  due 
to  internal  dissipation.  It  also  means  that  the  propa¬ 


gating  flexural  waves  reflect  off  the  ends  of  the 
numerical  ice  beam,  with  no  loss  to  the  fluid  or 
atmosphere.  One  of  these  reflected  waves  is  evi¬ 
dent  in  Figure  7a  at  about  time  0.25  s. 

Despite  the  shortcomings  noted  above,  the  nu¬ 
merical  model  shows  go<^  agreement  with  the 
experimental  data.  Furthermore,  for  impulsive 
force  measurements,  only  the  ice  response  during 
the  impact  is  of  concern.  This  is  the  portion  of  the 
event,  temporally,  where  the  numerical  model 
shows  the  best  agreement  with  the  experimental 
data. 


EXPERIMENTAL  RESULTS 

A  typical  impact  force/time  record  is  shown  in 
Figure  10.  The  X-direction  (along  the  axis  of  im¬ 
pact)  force  is  on  the  top,  the  Y-direction  (transverse 
to  impact)  is  in  the  middle  and  the  Z-direction 
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a.  At  the  impact  end. 


Figure  8.  Comparison  of  numerical  and  experimental  beam  displacement  histories  for  an  8-m-long  beam. 


(vertical)  is  the  bottom  curve. 

An  upward  force  on  the  bridge  is  represented 
by  a  downward  Z-force  trace.  The  bottom  line 
represents  the  carriage  speed.  It  was  included  to 
ccmfirm  duit  the  c2UTiage  speed  was  essentially 
constant  before  die  impact  event. 

AlookatFigure  10  is  informative  to  understand 
the  impact  event.  Hie  X-force  trace  (top)  shows, 
just  after  the  impact  peak,  Uiat  there  is  a  small  dip 
in  the  trace  where  the  force  falls  below  its 
pre-impact  value,  as  the  beam  reboimds  off  the 
face.  This  dip  may  be  caused  by  bridge's  rebound 
after  die  impact,  or  because  there  is  a  momentary 
pause  before  the  water  fills  in  behind  the  rebound¬ 
ing  beam.  The  trace  then  flattens  as  the  bridge 
continues  through  the  water,  until  catching  up  to 
the  rebounded  ^am.  Only  the  peak  values  of  the 
first  impact  are  used  in  this  report.  The  peak  value 
was  measured  from  an  average  value  calculated 
several  seconds  before  the  impact. 


The  Z-fbrce  trace  was  similar  with  the  exception 
of  the  model's  behavior  between  die  first  and 
second  impacts.  It  can  be  seen  from  die  Z-force 
trace  there  is  a  sinusoidal  vertical  forcing  diat 
appears  after  the  first  impact.  This  may  be  caused 
by  the  surging  of  die  water  as  it  fills  the  void  left  by 
die  beam  after  the  first  impact  or  perhaps  it  is  die 
free  response  of  the  bridge  after  the  impact 

The  Y-force  is  near  zero  as  it  should  be  in  a 
properly  aligned  impact.  Some  vibration  in  this 
direction  is  evident  shordy  after  the  X  direction 
peak  impact. 

The  data  on  maximum  total  ice  force  for  difid> 
ent  beam  thicknesses  and  lengths  are  presented 
graphically  in  Figures  11-14,  beginning  widi  die 
measurements  at  the  lowest  impact  velocities.  The 
magnitude  of  the  total  force  was  calculated  by 
adding  die  horizontal  and  vertical  compcments 
vectoiially.  The  tabulated  values  are  given  in  Ap¬ 
pendix  B. 
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Equation  9  was  used  to  obtain  prototype  values 
of  impact  force  from  the  results  of  the  model  tests 
and  these  values  are  also  given  in  Appendix  B. 

Figures  11-14  show  that  the  peak  impact  force 
increases  monotonically  with  increasing  velocity 
and  beam  length,  imtil  finally  exhibiting  a  ten¬ 
dency  to  level  off  and  approach  a  constant  value 
for  long  beams  and  higher  velocities.  This  phe¬ 
nomenon  is  quite  evident  from  the  test  results 
using  thin  beams.  As  beam  thickness  increases  the 
leveling  off  effect  is  not  significant  until  the  length 
or  velocity  increases  to  a  threshold  level.  An  in¬ 
spection  of  the  X  and  Z  forces  separately  indicates 
that  this  trend  holds  in  both  those  directions  as 
well  as  for  the  total  force. 

A  long  thin  beam  would  exhibit  this  trend  at  a 
lower  velocity  than  a  short  thick  beam  because  of 
a  thin  beam's  tendency  to  absorb  energy  through 
flexure,  while  a  short  thick  beam  acts  more  as  a 
rigid  body,  which  is  consistent  with  the  data.  The 


results  from  the  thickest  beams  did  not  exhibit  this 
trend  at  a  velocity  of  0.1524  m/s  regardless  of  the 
beam  length.  Some  reduction  of  force  barely  mani¬ 
fests  itself  at  0.60%  m/s  and  can  be  seen  rather 
well  at  1.0668  m/s.  Unfortunately  data  for  a  2.1-m- 
long  beam  at  1.524  m/ s  are  not  available,  and  this 
makes  it  difficult  to  see  the  trend  in  Figure  14. 

In  the  absence  of  friction,  bridge  inertia,  and 
any  effect  from  fluid  drag,  the  X  and  Z  component 
forces  would  be  equal  in  magiutude.  In  the  experi¬ 
ments,  however,  ^e  X-forces  are  on  the  average 
39%  greater  than  the  Z-forces.  Figure  15  is  a  graph 
of  X-force  versus  Z-force.  The  solid  line  is  the  best 
fit  regression  line  through  the  points.  If  the  differ¬ 
ence  in  the  component  forces  is  due  to  friction,  the 
coefficient  of  friction  (p)  can  be  calculated  by 

p  =  tan  ( (1))  (14) 

where  ^  is  the  friction  angle.  For  the  data  (}» is  8.30, 
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FigureW.  Typicalplotofforcex?s.  time  datafrom  the  experiments. 
In  this  case,  the  initial  impact  occurs  at  0.64  s. 
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Figure  13.  Average  peak  impact  force,  vs.  beam  length  and  beam  thickness,  for  the  replicated 
experiments,  at  1.0668  m/s. 


Figure  14.  Average  peak  impact  force,  vs.  beam  length  and  beam  thickness,  for  the  replicated 
experiments,  at  1.524  m/s. 
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Figure  15.  Maximum  values  ofF^vs.  F^,  showing  a  constant  ratio  ofF^^  be  approximately 
1.34.  This  yields  a  friction  angle  of  8.3°  and  a  Action  co^cient  f=FyF„)  of  0.15. 


which  gives  a  value  of  |x  of  0.15.  This  value  of  |xhas 
been  shown  to  be  possible  between  ice  and  stain¬ 
less  steel,  especially  when  the  stainless  steel  rough¬ 
ness  was  increased  (Forland  1985).  In  the  same 
report  it  was  also  shown  that  p  could  increase  with 
velocity. 

The  trend  inthedataisnot  uniform  and  exhibits 
a  rough  correlation  with  velocity.  For  any  given 
beam  size  the  X-force  is  closest  to  the  Z-force  at 
0.1524  m/s.  In  some  instances  the  Z-force  is  actu¬ 
ally  greater  than  the  X-force  at  this  velocity. 

In  general,  as  the  velocity  increases  the  X-force 
becomes  increasingly  larger  than  the  Z-force.  The 
notable  departures  from  this  trend  are  the  beams 
with  lower  length  to  thickness  ratios.  During  the 
impacts  described  earlier,  it  was  easier  for  short, 
thick  beams  to  rotate,  in  comparison  to  longer 
beams.  The  resulting  movement  in  the  Z  direction 
required  the  beam  to  displace  the  water  beneath  it. 
At  impact  these  unstable  beams  reacted  along  the 
path  of  least  resistance  (i.e.,  in  the  X  direction)  and 
therefore  generated  less  Z-force  against  the  45° 
sloped  surface.  This  is  discussed  later  in  this  re¬ 
port. 

It  was  speculated  that  the  increase  in  X-force 
might  be  associated  with  the  increase  in  viscous 


drag  caused  by  the  sudden  acceleration  of  the 
beam  through  the  water  at  impact.  To  obtain  the 
drag  force  associated  with  the  sudden  acceleration 
of  a  flat  plate,  the  shear  force  per  unit  area  is  first 
calculated  by 

^zx  =  -^^  (15) 

where  =  the  shear  stress  in  the  X  direction 
caused  by  fluid  in  the  Z  direction 
(perpendicular  to  the  plate) 
p  =  the  fluid  viscosity 
Vx  =  the  X-component  of  the  fluid  velocity 
(Bird  et  al.  1960). 

In  order  to  do  this  calculation,  an  estimate  of  the 
thickness  of  the  fluid  momentum  boundary  layer 
is  needed.  This  boundary-layer  thickness  5  is  ap¬ 
proximated  by 

6  =  4  Vw  (16) 

where  t  is  time  and  v  is  the  fluid  kinematic  viscos¬ 
ity 

u  =  — 
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Table  1.  Fluid  drag  forces  associated  with  the  sudden  acceleration  of  a 
flat  plate  for  each  size  ice  beam  calculated  at  the  highest  and  lowest 
experimental  velocities.  The  last  column  is  the  percentage  of  average  X- 
force  the  drag  force  comprises. 


Beam 

length 

(m) 

Velocity 

(m/s) 

Beam 

thickness 

(m) 

Fluid 

drag 

(N) 

Percentage 
of  average 
X-Jorce 

3.05 

0.1524 

0.0254 

11.37 

4 

3.05 

1.524 

0.0254 

113.8 

14 

3.05 

0.1524 

0.1524 

11.37 

1 

3.05 

1.524 

0.1524 

113.8 

1 

2.134 

0.1524 

0.0254 

7.9S 

3 

2.134 

1.524 

0.0254 

79.61 

10 

2.134 

0.1524 

0.1524 

7.95 

1 

2.134 

1.524 

0.1524 

79.61 
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with  p  the  fluid  density  (Bird  et  al.  1960).  It  is 
evident  from  eq  16  that  the  boundary  layer  thick¬ 
ness  is  proportional  to  the  square  root  of  time.  In 
the  calculations  of  the  fluid  drag  during  the  impact 
events,  the  rise  time  of  the  force  was  used  in  this 
calculation.  The  rise  time  is  the  time  from  the 
begiruiing  of  the  impact  to  the  time  of  maximum 
peak  force  measured  from  the  experimental  force 
time  histories.  This  time  was  relatively  consistent 
and  an  average  of  0.05  s  was  used  in  eq  16.  This 
gave  a  botmdary  layer  thickness  of  0.01  m.  To 
parametrically  bound  the  values,  the  shear  stress 
was  calculated  for  the  highest  velocity  of  1.5  m/s 
and  lowest  velocity  of  0.15  m/s.  The  stress  values 
from  these  calculations  were  then  multiplied  by 
the  beam  surface  area  to  arrive  at  a  final  drag  force 
value  for  a  given  size  beam  at  each  of  the  above 
velocities.  TTiis  value  was  then  divided  by  the 
average  X-force  value  at  the  same  conditions  to 
determine  the  percentage  of  the  X-force  that  could 
possibly  come  from  drag.  The  results  are  listed  in 
Table  1. 

Table  1  shows  that  this  drag  force  accounts  for 
less  than  15%  of  the  total  X-direction  force  in  the 
best  case.  The  surface  roughness  of  the  ice  is  not 
considered  in  this  calculation.  An  ice  beam  is 
probably  rougher  and  might  generate  greater  shear 
stress  than  what  the  flat  plate  model  provides,  but 
it  is  unlikely  that  it  could  account  for  the  difference 


between  15%  and  the  measured  39%  difference. 

During  the  experiments  it  was  noticed  that  the 
bridge  face  had  b^ome  dented.  This  denting  could 
have  changed  the  local  geometry  enough  to  present 
a  more  vertical  face  to  the  impacting  ice  thus 
causing  higher  X-direction  forces. 

Another  possible  explanation  for  the  higher  X- 
forces  might  be  a  dynamic  interaction  of  the  model 
during  beam  impact.  The  configuration  of  the 
model  and  transducer  relative  to  where  the  ice 
impacted  the  model  could  cause  this  interaction  to 
show  up  more  in  the  X-direction  than  in  the 
Z-direction. 


RESULTS  AND  DISCUSSION 

Figure  16  shows  a  typical  force  versus  time 
record  for  a  1.219-m-long,  152-iiun-thick,  1.37-m- 
wide  beam  from  the  numerical  program.  Since  the 
bridge  is  modeled  as  a  massless  structure,  the  force 
versus  time  record  is  symmetric  with  no  effect  of 
bridge  inertia  evident  in  the  impact.  The  force 
values  used  in  the  comparisons  below  were  the 
maximum  impact  force  corresponding  to  the  apex 
of  the  trace. 

Figure  17  shows  comparisons  between  the  nu¬ 
merical  peak  impact  force  and  the  average  experi¬ 
mental  peak  impact  force,  versus  ice  beam  length. 
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Figure  16.  A  typical  numerical  force  vs.  time  plot  from  an  impact  of  a 
1.219-m-long  beam  at  0.61  m/s. 
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Figure  1 7.  Comparisons  of  average,  experimental,  peak  impact  force,  and  numerical,  peak  impact  force  vs.  ice 
beam  length,  at  constant  velocities,  for  a  given  thickness  ice  beam. 
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Figure  18.  Comparisons  of  average  experimental  peak  impact  force,  and  numerical  peak  impact  force  vs.  impact 
velocity,  at  constant  beam  lengths,  for  a  given  thickness  ice  beam. 


at  constant  velocities,  for  each  experimental  ice  In  both  figures,  the  numerical  force  is  higher 

thickness.  It  can  be  seen  here  that  the  numerical  than  the  experimental  force  for  the  thinner  ice  at  all 

and  experimental  forces  exhibit  the  same  tendency  but  the  lowest  velocities,  becoming  closer  to  the 

to  approach  a  constant  value  as  beam  length  and  experimental  force  as  the  ice  thickness  increases, 
velocity  increase.  This  is  understandable,  since  numerically  the  ice  is 

In  Figure  18,  comparisons  between  numericcd  treated  as  an  elastic  body,  with  no  failure  criteria, 
peak  impact  force  and  average  experimental  peak  widi  the  energy  completely  conservative.  In  the 

impact  force,  versus  impact  velocity  at  constant  experiments,  the  thiimer  ice  was  easier  to  damage, 

beam  length,  are  shown  for  each  experimental  ice  especially  as  the  velocities  increased,  which  means 

thickness.  Here  it  appears  that  the  impact  force  theexperimentalforceswouldnotbeashighasthe 
from  both  the  numerical  and  experimental  0.305-  numerically  modeled  forces.  Most  of  the  damage 
m-long  beams,  when  plotted  against  velocity,  ex-  in  our  experiments  took  the  form  of  either  cracks 

Mbits  a  different  trend  than  that  of  the  other  beam  propagating  back  from  the  nose  of  the  beam,  or 

lengths.  Generally,  the  slope  of  the  curve  for  the  pieces  of  the  beam  fracturing  and  breaking  off. 

0.305-m-longbeam  is  flatter  than  those  of  the  other  Appendix  B  lists  the  damage  noticed  in  the  experi- 

beam  lengths,  and  the  force  magnitudes  appear  to  ments. 

be  noticeably  lower.  The  forces  generated  during  the  ice  impact 
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Figure  19.  During  impact,  only  the  initial  portion  of  an  ice  beam  will  contribute  to  the  impact force 
as  illustrated  above.  The  size  (f  the  portion  that  contributes  to  the  impact  will  depend  upon  the  length 
and  characteristic  length  of  the  ice  beam. 


events  depend  upon  factors  such  as  the  mass  of  the 
bodies,  the  impact  velocity,  the  post-impact  ve¬ 
locities  of  the  bodies,  the  internal  dissipation  of 
energy  in  the  bodies,  and  the  energy  lost  to  the 
supporting  fluid  during  the  impact. 

For  the  experimental  impacts  where  the  ice 
beams  broke  and  fractured,  the  impulsive  force 
recorded  would  be  lower  than  what  a  pure  elastic 
body  would  have  generated  because  of  the  energy 
used  during  fracturing.  In  the  impacts  where  the 
beams  did  not  break,  a  linear  force  increase  ac¬ 
cording  to  Newton's  second  law  was  expected.  As 
Figure  17  shows,  however,  the  forces  tended  to 
approach  a  constant  value  with  increasing  beam 
length.  This  trend  is  also  evident  in  the  numerical 
data,  so  it  is  not  solely  related  to  energy  used  in 
deunaging  the  ice. 

This  trend  could  be  related  to  the  characteristic 
length  of  the  ice,  which  would  imply  that  impact 
force  will  approach  a  constant  value  when  the 
beam  length  is  more  than  some  multiple  of  the 
characteristic  length.  Any  length  increase  beyond 
this  multiple  will  not  contribute  to  the  impact 
force,  as  only  the  front  part  of  the  beam  will  react 
during  the  impact.  This  is  illustrated  in  Figure  19. 
If  this  supposition  is  correct  then  it  should  be 
possible  to  nondimensionalize  and  sc2de  the  data 
with  respect  to  the  characteristic  length. 

The  characteristic  lengths  of  the  ice  beams  ate 
listed  in  Table  2.  These  were  calculated  using  eq  12 
and  13. 

Figure  20  shows  comparisons  between 
nondimensionalized  numericsd  and  experimental 
impact  force,  versus  beam  length  divided  by  char¬ 
acteristic  length  (lyif)  at  constant  velocities.  The 
force  is  nondimensionalized  by  dividing  by  the 
effective  volume  of  the  beam,  which  is  calculated 


using  the  beam's  characteristic  lengdt,  multiplied 
by  the  acceleration  of  gravity  and  the  demity  of 
water.  This  results  in  a  nondimensional  force  v^ue 
that  has  the  effective  mass  of  the  beam  incorpo¬ 
rated  in  it.  The  effect  of  beam  duckness  is  also 
illustrated  within  each  figure. 

The  nondimensioiud  numerical  force  does  tend 
to  approach  a  constant  value  with  increasing  LA^ 
as  Figure  20  shows.  The  forces  from  the  0.305-m- 
longbeams  (low  IVI^  values)  ate  consistenfiy  higher 
dum  the  odter  length  beams  at  similar  velocities. 
The  experimental  forces  overall  are  lower  than 
what  is  predicted  by  the  munerical  data.  This  is 
expected  because  theexperimentalice  beams  could 
be  damaged  during  impact  while  the  numerical 
beams  are  modeled  as  pure  elastic  bodies.  The 
difierences  between  the  numerical  predictions  and 
measured  experimental  forces  were  less  at  die 
lowest  velocities,  where  damage  to  the  experi¬ 
mental  beams  was  less  likely  to  occur. 

A  comparison  of  the  experimental  and  muneri¬ 
cal  nondimensioi\al  force  data  shows  overall  a 
much  flatter  trend  of  force  vs.  L/L^  in  the  experi¬ 
mental  data.  The  force  increases  more  steeply  at 
lower  L/Lc  values  and  then  tends  to  flatten  out  with 
increasing  L/L^,  in  the  experimental  data,  for  all 
but  the  25-mm-thick  ice.  The  force  then  tends  to 
approach  a  constant  V2due  L .  much  the  same  man¬ 
ner  as  the  munerical  force,  except  with  a  lower 


Table  2.  Beam  thickness  and  associated  character¬ 
istic  length  for  the  test  beams. 


Thickness  (m)  0.0254  0.0508  0.1016  0.1524 

Lc(in)  0.70  1.18  1.98  2.68 
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Figure  20.  Comparisons  ofnondimensional  experimental  and  numerical  impact  force,  vs.  beam  length  divided  by 
characteristic  length  (L/L^J,  at  constant  beam  thickness,  for  a  given  impact  velocity.  The  length  divided  by 
thickness  ratio  (L/h)  of  the  ice  beam  is  beside  the first few  experimental  points  that  do  not  follow  the  trend  of  high 
nondimensional  force  at  low  L/L^  ratios. 


magnitude.  The  numerical  data  also  follow  this 
trend  where  they  initially  increase  relatively 
quickly  at  low  L/L<.  ratios,  then  appear  to  asymp¬ 
totically  approach  a  constant  force. 

The  above  comparison  suggests  that  the  short¬ 
est  beams  are  behaving  differently  than  the  longer 
beams  during  impact.  The  data  show  a  more  steeply 
rising  nondimensional  force  for  the  short  beams, 
which  then  drops  off  with  increasing  L/L^  ratios. 
This  difference  is  a  function  of  the  beeun  length/ 
thickness  ratio  (L/h),  which  determines  how  the 
beam  reacts  during  impact.  A  very  small  length/ 
thickness  ratio  would  mean  that,  in  profile,  the 
beam  would  approach  a  square  shape,  which 
would  be  unstable  and  easy  to  roll  over  at  impact. 
The  shortest  ice  beams  would  then  roll  vertically 
over  their  upstream  bottom  comer  as  they  are 
shoved  backwards  into  the  fluid  behind  the  beam 
during  the  impact.  This  behavior  would  not  be  as 


likely  in  the  numerical  simulation,  because  there  is 
no  friction  between  the  bridge  face  and  ice,  and 
therefore  no  mechanism  to  push  the  beam  back¬ 
wards.  Furthermore,  there  is  also  no  "fluid"  to 
hold  the  upstream  bottom  comer  when  the  beam 
is  pushed  backwards.  In  the  numerical  model, 
these  short  beams  might  behave  like  a  rigid  body 
during  impact,  which  would  explain  the  relatively 
higher  peak  force  in  the  shortest  beams,  vs.  the 
impacts  against  longer  or  thinner  beams  where 
some  of  the  energy  went  into  flexing  the  ice.  This 
would  be  similar  to  the  energy  required  to  com¬ 
press  a  spring.  This  scenario  would  not  be  as  likely 
in  the  short,  thin  ice  beams,  since  the  L/h  ratio  of  a 
short,  thin  beam  is  not  as  large  as  that  of  the  thick, 
short  beams  and  thus  not  as  prone  to  rolling  dur¬ 
ing  impact. 

The  trends  in  the  numerical  data  support  the 
argument  that  as  the  L/L^.  ratio  of  a  beam  becomes 


20 


Figure  21.  Nondtmensioml  experimental  and  numerical  force  vs.  nondimensional 
velocity.  The  arrow  at  approximately  1.1  on  the  velocity  scale  indicates  where  the 
majority  of  the  experimental  beams  started  to  be  visibly  damaged  at  impact. 


larger  the  force  will  approach  a  constant  value. 
The  impact  velocity  will  affect  the  magnitude  of 
this  value.  The  experimental  data  also  tend  to 
approach  a  constant  value  with  increasing  UL^ 
ratios,  but  it  is  lower  overall  than  what  the  numeri¬ 
cal  data  suggest.  This  could  be  because  as  the 
velocity  increases,  the  experimental  beams  are 
being  damaged,  and  thus  cannot  sustain  as  high 
an  impulsive  force  as  the  numerical,  purely  elastic 
beam.  This  is  illustrated  in  Figure  21,  which  plots 
nondimensional  force  vs.  velocity,  for  both  the 
numerical  and  experimental  data.  In  the  figure, 
the  points  representmg  the  0.305-m  beams  are 
separate  from  the  other  data.  The  upward  point¬ 
ing  arrow  on  the  velocity  axis  indicates  the  lowest 
velocity  at  which  the  majority  of  the  experimental 
beams  became  damaged.  In  the  figure,  ^e  numeri¬ 
cal  and  experimental  data  start  to  diverge  close  to 
the  point  where  the  experimental  beams  start  to 
break.  As  the  velocity  increases  the  discrepancy 
becomes  greater,  while  at  the  lower  velocities,  the 
numerical  and  experimental  data  compare  well. 

In  Figure  21  the  nondimensional  forces  from 
the  experimental 0.305-m  beams  are  below  the  rest 
of  the  experimental  data,  and  the  forces  from  the 
numerical  0.305-m  beams  are  also  below  the  other 
numerical  data.  This  trend  from  the  0.305-m  long 


beams  is  evident  at  nearly  every  impact  velocity 
and  again  supports  the  suggestion  that  the  0.30^ 
m  beams  are  behaving  differently  than  the  other 
beams.  The  other  differences  between  the  numeri¬ 
cal  and  experimental  data  for  beams  of  all  lengths 
could  be  explained  by  the  fact  that  the  experimen¬ 
tal  beams  are  breaking  while  the  numerical  beams 
are  modeled  as  perfectly  elastic  bodies  with  infi¬ 
nite  breaking  strength. 

Numerical  experiments  were  performed  to  in¬ 
vestigate  the  effect  of  changing  material  proper¬ 
ties  between  the  bridge  model  and  ice  beam.  Fig¬ 
ure  22a  shows  the  minimal  effect  that  the  change  in 
ice  elastic  modulus  had  upon  impact  force  for  the 
152-mm-thick  ice.  A  much  greater  effect  can  be 
seen  in  Figure  22b  where  the  effect  of  changing  the 
structure  stiffness  is  illustrated.  At  the  time  of  the 
experimental  study,  modeling  the  structure  as  a 
rigid  body  was  thought  to  be  acceptable,  and  that 
this  approximation  would  give  reasonable,  if  not 
somewhat  high,  impact  values.  Figure  22b,  how¬ 
ever,  illustrates  the  importance  of  knowing  the 
structure  stiffness  during  ice  impact  events.  It  is 
evident  that  a  change  in  structure  stifhiess  would 
greatly  affect  the  impact  force  values  and  an  as¬ 
sumption  that  the  structure  is  rigid  could  give 
unacceptable  predictions. 
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a.  Plot  showing  the  effect  of  changing  the  ice  elastic  modulus from  2.0 
tolO.OGPa. 
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b.  Plot  showing  the  ^ect  of  changing  the  structure  stiffness,  from  250  x 
lO^toSOOxlO^N/m. 

Figure  22.  Plots  of  force  vs.  xfelocity,  at  constant  beam  lengths,  from  the 
numerical  simulation.  These  compare  the  effect,  on  impact  force,  of 
changing  the  ice  elastic  modulus  vs.  changing  the  structure  stiffness. 


PROTOTYPE  FORCES 

Some  of  the  prototype  force  values  that  were 
determined  from  this  study  are  disturbing  in  that 
they  seem  to  be  extremely  high.  However,  it  should 
be  remembered  that  the  model  ice  beam  sizes 


represent  a  range  of  prototype  sizes  from  1.4  m  to 
79  m.  These  beams  were  propelled  at  a  top  proto¬ 
type  speed  of  3.4  m/s.  The  larger  beams  represent 
a  substantial  momentum  change  and  thus  high 
forces  at  the  highest  velocity  tested. 

To  check  the  findings,  a  comparison  of  the 
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results  of  this  study  with  those  obtained  in  the 
field  would  be  necessary.  For  the  field  data  to  be 
useful  to  us,  the  data  needed  to  include  stream 
velocities,  peak  impact  forces,  impact  times  and 
the  size  of  die  beam  that  produced  the  force.  Some 
studies  of  ice  booms  that  were  investigated  re¬ 
corded  the  first  two  pieces  of  information  (Perham 
1977, 1978, 1983)  but  not  the  impact  times  nor  the 
size  of  the  ice  floes  creating  the  force.  Some  of  the 
peak  forces  recorded  by  Perham  (1977)  were  in  the 
150-  to  730-kN  range. 

However,  it  is  difficult  to  compare  the  ice  boom 
force  measurements  with  the  impact  data.  The 
measured  peak  force  from  any  impact  is  depen¬ 
dent  upon  the  dmation  of  the  impulse.  An  ice 
boom  can  absorb  energy  through  deflection  of  the 
boom  and  therefore  may  show  a  much  lower  peak 
force  for  any  given  impulse  value.  The  peak  force 
in  the  modeled  situation  would  probably  be  higher 
for  any  similar  floe  that  strikes  an  ice  boom,  since 
the  bridge  model  is  more  rigid  than  a  boom.  Fur¬ 
thermore,  many  of  these  ice  boom  forces  tended  to 
be  caused  by  wind  and  water  drag  on  huge  ice 
sheets  as  opposed  to  discrete  floes  striking  the 
booms. 

A  report  on  ice  force  measurements  at  a  bridge 
pier  (Sodhi  and  Gagnon  1989,  Sodhi  et  al.  1983), 
however,  gives  force  measiuements  attributed  to 
relatively  discrete  ice  floes  impacting  an  instru¬ 
mented,  vertical  sided,  bridge  pier.  By  direct  ob¬ 
servation,  the  river  velocity  was  estimated  to  be  1.2 
to  1.5  m/s,  with  some  of  the  larger  floes  that  were 
observed  during  impact  estimated  to  be  19  to  25 
m^  large.  The  forces  recorded  were  predominantly 
in  the  25-  to  90-kN  range  with  many  above  100  kN. 
The  lower  forces  presmnably  came  from  the  smaller 
brash  ice  pieces  striking  the  pier.  Several  impacts 
were  from  approximately  130  to 200  kN.  The  high¬ 
est  actual  force  recorded  was  310  kN. 

Our  study  had  three  model  beam  sizes  that 
were  close  to  this  protot)rpe  range:  17.8, 21.7  and 
25.5  m2.  Our  modeled  prototype  velocity  closest  to 
the  above  velocity  range  is  1.4  m/s.  Our  modeled 
protot3rpe  forces  for  these  data  range  from  137  to 
306  kN,  which  gives  good  agreement  with  the 
larger  floe  forces  measured  on  vertical,  rigid  bridge 
piers.  The  prototype  forces  for  all  of  the  data  are 
listed  in  Ap{>endix  B. 


CONCLUSIONS 

A  floating  besun  was  subjected  to  vertical  im¬ 
pulsive  loads  and  the  resultant  beam  displace¬ 


ments  were  measured.  The  displacements  mea¬ 
sured  in  the  experiments  were  modeled  using  a 
finite  element  analysis  that  designated  the  ice  eis  a 
linear  elastic  material  The  dynamic  effect  of  the 
fluid  foimdation  was  calculate  in  a  separate  finite 
element  program  and  added  to  the  numerical 
beam  in  the  same  maimer  as  an  added  mass  ma¬ 
trix.  Good  agreement  between  the  numerical  and 
experimental  results  confirmed  the  validity  of  the 
el^tic  designation  and  the  fluid  added  mass. 

A  series  of  experiments  were  performed  to 
assess  the  impact  force  generated  by  freshwater 
ice  beams  of  four  different  lengffis  and  thicknesses 
at  four  velocities.  The  impact  forces  from  the  ex¬ 
periments  were  munerically  modeled  using  finite 
element  analysis,  with  an  elastic  material  designa¬ 
tion  for  the  ice  and  die  fluid-added  mass  devel¬ 
oped  earlier.  The  conclusions  ffom  the  work  are: 

1.  The  finite  element  method  used  in  modeling 
ice/structure  impacts  allows  one  to  easily  model 
complex  geometries,  vaiyingimpactvelocities  and 
different  material  properties.  In  the  model,  the  ice 
is  characterized  as  an  elastic  body,  with  a  nonlin¬ 
ear  foimdation  modulus  to  account  for  the  static 
fluid  forces,  and  a  fluid  influence  coefficient  ma¬ 
trix  to  accoimt  for  the  hydrodynamic  fluid  forces. 
The  impact  forces  obtained  from  the  model  gener¬ 
ally  compare  well  with  the  experimental  data  for 
the  thickest  ice  beams  and  at  the  lower  velocities 
for  the  other  thicknesses.  The  largest  discrepancies 
appear  in  the  forces  associated  with  the  very  short 
(0.305-m  long)  beams.  There  are  also  differences 
evident  when  the  ice  fails,  which  are  not  accounted 
for  in  oiur  model 

2.  A  failure  criteria  needs  to  be  developed  for 
the  numerical  model  to  properly  predict  the  im¬ 
pact  force  when  the  ice  beams  break. 

3.  The  impact  forces  caused  by  the  short  beams 
are  not  modeled  well  in  our  present  numerical 
model  As  discussed  in  this  report,  the  inclusion  of 
fluid  drag  against  the  bottom  and  upstream  face  of 
the  ice  beams  along  wifii  the  addition  of  friction  at 
the  bridge/ice  interface  could  correct  the  short 
beam's  vertical  rotation  after  impact  and  possibly 
bring  the  numerical  predictions  closer  to  the  ex¬ 
perimental  results. 

4.  Ice  beam  impact  forces  appear  to  approach  a 
constant  value  with  increasing  beam  length  for  a 
given  impact  velocity.  This  is  thought  to  be  a  result 
of  the  energy  absorbed  in  the  flexing  of  the  beam 
under  an  impact  load,  and  is  described  by  using 
the  ice's  characteristic  length.  The  experimental 
results  appear  to  agree  with  the  numerical  data  up 
until  the  ice  is  damaged.  The  thinnest  ice  beams 
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exhibit  the  same  relative  trends  as  the  numerical 
data  at  all  velocities.  It  should  be  possible  to  arrive 
at  a  theoretical  maximum  force,  at  a  given  velocity, 
for  long  slender  beams,  based  upon  a  ratio  of 
length  divided  by  characteristic  length.  In  the  data, 
a  length/characteristic  length  ratio  of  3  to  5  is 
where  the  forces  appear  to  approach  a  constant 
value,  but  more  data  for  longer  beams  are  needed 
to  confirm  this  upper  limit. 

5.  The  structure  stiffness  is  important  in  deter¬ 
mining  the  impact  force  in  ice/ structure  interac¬ 
tions.  The  shape  and  magnitude  of  the  impulsive 
force  depends  heavily  upon  the  reaction  of  the 
structure  to  the  impact. 

6.  The  prototype  forces  calculated  from  the 
experiments  appear  to  be  in  agreement  with  lim¬ 
ited  available  data  of  actual  force  measurements 
for  those  cases  where  the  size  and  velocities  of  the 
impacting  floes  were  known. 
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APPENDIX  A:  PRELIMINARY  RESULTS  AND  ANALYSIS 


A  typical  preliminary  impact  force/time  record  is  shown  in  Figure  Al.  The  X 
direction  is  the  top  trace,  the  Y  direction  is  the  middle  trace  and  the  Z  direction  is  the 
bottom  trace.  A  downward  Z-force  trace  represents  an  upward  force  on  the  bridge. 
The  impact  occurs  at  about  time  1  second  and  it  is  evident  that  the  X-force  is  larger 
than  the  Z-force.  The  face  of  the  bridge  is  at  a  45®  angle  to  the  beam  which  should 
cause  nearly  equal  X  and  Z  resultant  forces  on  the  bridge.  In  addition,  there  is  also 
a  substantial,  oscillating,  tratrsverse  Y-force.  If  the  beam  is  lined  up  perpendicular 
to  the  bridge  face  and  strikes  it  squarely  then  the  amoimt  of  force  transmitted 
transverse  to  the  face  (Y-direction)  will  be  minimal,  probably  mostly  due  to 
vibration  in  the  model. 

Since  it  is  difficult  to  align  a  wide  beam  to  strike  the  bridge  squarely,  an  off-center 
impact  would  provide  a  larger  moment  arm  to  rotate  the  beam  about  the  vertical 
axis.  Figure  A2  shows  various  configurations  of  beam  shapes  that  were  used  to 
determine  if  the  scatter  in  the  data  was  due  to  the  beams  impacting  the  bridge 
off-center. 

Figure  A3  is  a  graph  of  X  and  Z  impact  forces,  plotted  against  the  relative  mass 
of  each  beam.  The  largest  beam  had  a  relative  mass  equal  to  one.  The  beam  mass  was 
changed  by  changing  the  width  of  the  beams  and  by  removing  material  for  the 
different  nose  shapes. 

If  the  impacts  are  elastic  cmd  no  energy  is  being  lost  to  beam  rotation  or  breakage, 
then  the  forces  recorded  should  follow  a  linear  trend  against  mass,  and  pass  through 
the  origin.  In  addition,  if  the  experiments  are  consistent,  the  scatter  of  the  points  at 
any  relative  mass  should  be  minimal.  It  can  be  seen  that  the  1 .37-m- wide  beams  with 
a  full  width  nose  had  a  wide  scatter  in  both  the  X  and  Z-forces.  The  next  largest 
scatter  came  from  the  GMC  beams  (see  Fig.  A2).  These  beams  had  a  0.76-m-wide 
nose  that  angled  back  to  the  main  body  of  the  beam  at  a  45°  angle  (Fig.  A2).  This  nose 
is  still  relatively  wide,  which  probably  contributed  to  much  of  the  scatter.  In 
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Figure  Al.  Typical  force  vs.  time  record  from  the  preliminary  experi¬ 
ments. 
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Figure  A3.  Plot  cfP^  relative  mass. 


addition,  the  Z  direction  force  imparts  a  bending  stress  to  the  beams,  causing  some 
to  break  in  tension,  at  the  top  surface.  The  long  and  relatively  narrow  beam  area 
behind  the  nose  would  present  less  material  to  absorb  the  biding  stresses  than 
would  a  full  width  beam.  This  might  have  contributed  to  the  wider  scatter  in  the 
forces  for  these  beams.  It  would  be  expected,  though,  that  if  the  beams  were 
breaking  prematurely  due  to  the  narrower  area  behind  the  nose,  the  average  force 
for  these  beams  would  be  lower  than  what  it  was.  It  is  evident  that  the  forces  from 
the  beams  with  the  narrower  noses  were  more  consistent. 

The  beams  labeled  OC  and  GOC,  in  Figure  A2,  were  modified  to  piuposely  hit 
off-center  and  rotate  during  impact.  Both  the  X  and  Z-forces  recorded  for  the  GOC 
beams  are  well  below  the  others.  These  beams  were  the  most  off-center  design  used, 
and  it  is  evident  that  a  lot  of  energy  was  expended  to  rotate  the  beams.  The  OC  beams 
did  not  contain  as  much  mass  to  one  side  as  the  GOC  beams  did  and  are  more 
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Figure  A4.  Plot  ofpj^  and  vs.  beam  relative  mass,  and  a  best  fit  linear  regression 
line,  through  the  data,  from  the  origin. 


indicative  of  a  full  width  beam  striking  slightly  out  of  alignment.  These  points  are 
below  the  majority  of  the  data,  with  the  exception  of  the  lower  scatter  from  the  full- 
width  beams.  This  is  consistent  with  the  speculation  that  the  lower  scatter  of  the  full- 
width  beams  is  caused  by  impact  with  a  slightly  skewed  beam. 

The  second  characteristic  looked  for  in  the  data  was  a  linear  trend  through  the 
origin.  It  can  be  seen  in  Figure  A3  that  beams  with  FW,  OC  and  GCXZ  shapes  fall 
below  any  clearly  defined  linear  trend  for  both  the  X-force  and  Z-force  data. 
Furthermore,  the  GMC  beams  were  excluded  because  of  the  scatter  and  potential 
breaking  problems  associated  with  the  longer  nose  area  of  these  beams. 

Figure  A4  shows  the  X  and  Z-force  1.07  MC,  0.76  MC,  1.37  MC  and  Unmod  HW 
data  that  are  retained,  with  the  best  fit  linear  regression  line  through  the  origin  that 
represents  each  set.  Table  A1  lists  the  model  and  correlation  coefficients  for  each. 

The  X-force  data  are  fairly  well  represented  by  a  linear  fit  model.  A  linear  fit 
through  the  origin  does  not  represent  the  Z-force  data  quite  as  well,  but  is  still 
reasonable. 

Table  A2  shows  the  sample  size,  average  force  and  standard  deviation  for  the 
component  forces  of  each  type  of  beam.  The  magnitude  of  the  standard  deviation 
compared  to  the  average  force  is  also  listed  and  can  be  used  as  an  indicator  of  the 


Table  Al.  Slope  and  correlation  coefficients  for  the 
regression  lines  shown  in  Figure  A4. 

Slope  Correlation  coefficient 

X-force  data  1653.85  0.9966 

Z-force  data  1275.91  0.9856 
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Table  A2.  Sample  size,  average  force,  standard  deviation,  and  standard  deviation  divided  by 
average  force,  of  the  component  forces,  for  each  type  of  beam  used  in  the  preliminary 
experiments. 

Floe  description 

HW 

HW 

FW 

FW 

OC 

OC 

Force  direction 

X 

Z 

X 

Z 

X 

z 

Sample  size 

26 

26 

21 

20 

4 

4 

Average 

1213.6 

799.8 

1238.5 

801.1 

1037.5 

642.7 

Standard  deviation 

420.4 

391.7 

470.3 

281.0 

124.0 

154.5 

%  of  avg  (S.DVavg) 

35% 

49% 

38% 

35% 

12% 

24% 

Floe  description 

GOC 

GOC 

1.07  MC 

1.07  MC 

.76  MC 

.76  MC 

Force  direction 

X 

Z 

X 

Z 

X 

Z 

Sample  size 

6 

6 

4 

4 

4 

4 

Average 

675.4 

329.2 

13%.7 

1076.4 

844.0 

673.9 

Standard  deviation 

58.6 

69.5 

70.3 

57.8 

41.9 

61.5 

%  of  Avg  (S.DVavg) 

12% 

21% 

5% 

5% 

5% 

9% 

Floe  description 

1.37  MC 

1.37  MC 

GMC 

GMC 

Force  direction 

X 

Z 

X 

Z 

Sample  size 

4 

4 

10 

10 

Average 

1623.5 

1382.2 

1541.2 

987.9 

Standard  deviation 

64.3 

72.7 

163.4 

162.8 

%  of  avg  (S.DVavg) 

4% 

5% 

11% 

17% 

repeatability  achieved  versus  beam  type.  The  modified  center  beams  (MC)  gave  the 
best  repeatability  as  is  evidenced  by  the  low  percentage  figures.  For  this  reason  the 
modified  center  beams  were  the  design  used  in  the  experiments. 

If  the  impact  strictly  follows  Newton's  2nd  law  relationship,  then  varying  the 
width  of  the  beams,  all  else  being  equal,  should  not  affect  the  impact  force,  beyond 
what  the  corresponding  change  in  mass  produces.  Figure  A5  shows  the  X  and  Z- 


Width  (m) 

Figure A5.  PlotofFj^andF^vs-  icebeamwidthfortheOJBMC,  1.07MCand 
1.37  MC  beams. 


Table  A3.  Slope  and  correlation  coefficients  for  the 
regression  lines  shown  in  Figure  A5. 

Slope  Correlation  coefficient 

X-foicedata  1211.34  0.9952 

Z-forcedata  988.37  0.9954 


forces  from  the  0.76  MC,  1.07MC  and  1.37MC  beams.  In  these  beams  only  the  width 
was  varied  while  the  nose  design  stayed  constant.  Although  the  data  did  not  show 
a  perfect  linear  force  increase  with  increasing  width,  it  appears  quite  good  as  is 
evidenced  by  the  be  the  best  fit  linear  regression  line  through  the  origin  for  each  set. 
Table  A3  shows  the  slopes  and  correlation  coefficients  for  the  data. 

In  use,  the  bridge's  displacement  will  change  depending  upon  the  load  it  is 
carrying.  It  was  not  expected  that  the  impact  force  would  be  greatly  affected  by  this 
change  but  it  was  prudent  to  check  it.  Figure  A6  is  a  graph  of  impact  force  vs  bridge 
displacement.  The  displacement  values  used  modeled  the  full-scale  range,  from  the 
no-load  displacement  up  to  the  maximum  expected  load  value.  A  higher  displace¬ 
ment  value  means  the  bridge  is  deeper  into  the  water.  There  does  not  appear  to  be 
a  significant  change  in  force  between  the  first  two  displacements.  However,  at  the 
maximum  displacement  the  X-force  appears  to  drop,  while  the  Z-force  rises. 

At  the  low  displacement  the  X-force  is  greater  than  the  Z-force  by  an  average 
138%.  It  also  appears  that  as  the  displacement  becomes  greater  the  two  component 
forces  approach  each  other  in  magnitude.  With  a  45®  sloped  face  one  would  expect 
the  components  to  be  close  to  equal  disregarding  the  effect  of  friction  or  a  dynamic 
interaction  of  the  model.  This  characteristic  of  higher  X-forces  is  noted  in  the  main 
data  and  discussed  in  the  report. 


Figure  A6.  PlotofFj^and  F^ds.  model  displacement.  The  beam  was  152 
mm  thick,  1.219  m  long  and  the  impact  velocity  was  0.6096  m/s. 
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APPENDIX  B:  TEST  DATA. 


Damage  codes  in  damage  column: 

O  -  No  visible  damage 
C  -  Beam  cracked 
B  -  Beam  broke 


Model  data 

Model  width  =  1.372  m 


Floe  Floe 

thick,  length  Velocity  X-force  Z-force  Total 
(m)  (m)  (m/s)  (N)  (N)  (N) 


0.1324 

ElIKIH 

916 

763 

1193 

i 

i 

01324 

1317 

1148 

1902 

• 

• 

0.1324 

1788 

1210 

2139 

• 

• 

01324 

1303 

1108 

1867 

• 

• 

0.6096 

7024 

4826 

8322 

• 

• 

0.6096 

6348 

4662 

8067 

• 

• 

0.6096 

4938 

3398 

3994 

• 

• 

1.0668 

8739 

6161 

10709 

• 

• 

1.0668 

8314 

6076 

10460 

• 

• 

1.0668 

8394 

8394 

11870 

• 

• 

1.0668 

9110 

6130 

10980 

• 

• 

1.324 

10920 

7744 

13388 

• 

T 

1.324 

11614 

8496 

14390 

• 

3.04S 

1.324 

10916 

8063 

13372 

• 

0.1324 

747 

614 

967 

m 

01324 

907 

712 

1133 

• 

• 

0.1324 

810 

810 

1149 

• 

• 

01324 

827 

827 

1170 

• 

• 

0.6096 

4393 

3087 

3371 

• 

• 

0.6096 

3107 

3278 

6068 

• 

• 

0.6096 

3022 

3136 

9921 

• 

• 

0.6096 

4493 

3149 

3487 

• 

T 

1.0668 

7638 

8962 

• 

2.134 

1.0668 

7228 

4466 

8497 

• 

1.219 

0.1324 

418 

387 

370 

• 

i 

0.1324 

334 

449 

698 

• 

• 

0.1324 

363 

472 

736 

• 

• 

0.6096 

2140 

1326 

2628 

• 

• 

0.6096 

2126 

1366 

2641 

• 

• 

0.6096 

2393 

1463 

2803 

• 

• 

0.6096 

2602 

1268 

2893 

• 

• 

1.0668 

4017 

2237 

4398 

• 

• 

1.0668 

4330 

2220 

4884 

• 

• 

1.0668 

4279 

2237 

4829 

• 

• 

1.0668 

4133 

2246 

4723 

• 

• 

1.324 

3182 

3309 

6149 

m 

1.324 

3760 

3492 

6736 

1.219 

1.324 

3843 

2918 

6333 
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Prototype  data 
Prototype  width  =  6.86  m 

Floe  Floe  Total 

thick.  length  Vel.  force 

Dam.  (m)  (m)  (m/s)  (kN) 


Damage  codes  in  damage  column: 


O  -  No  visible  damage 
C  -  Beam  cracked 
B  -  Beam  broke 

Floe  Floe 

thick,  length  Velocity  X-force  Z-force  Total 
(m)  (m)  (m/s)  (N)  (N)  (N) 


1  0.1524 

0.305 

0.1524 

89 

85 

123 

mm 

i 

0.1524 

83 

85 

120 

9 

9 

0.1524 

tnmm 

116 

164 

• 

E!12m 

tnmm 

120 

158 

9 

• 

0.6096 

476 

334 

581 

9 

• 

0.6096 

480 

400 

625 

• 

• 

0.6096 

431 

285 

517 

• 

• 

0.6096 

458 

280 

337 

• 

• 

1.0668 

916 

267 

954 

• 

9 

1.0668 

932 

307 

1000 

• 

9 

1.0668 

925 

569 

1086 

• 

• 

1.0668 

903 

436 

1003 

• 

9 

1.324 

1397 

472 

1474 

« 

9 

1.524 

1435 

276 

1480 

m 

T 

1.524 

1477 

440 

1541 

0.1324 

0.305 

1.324 

1708 

445 

1763 

0.1016 

3.04S 

796 

805 

1132 

mm 

i 

0.1524 

801 

850 

1167 

9 

9 

0.1524 

894 

930 

1290 

9 

9 

0.1514 

921 

952 

1324 

9 

9 

0.6096 

4097 

3545 

5418 

9 

9 

0.6096 

3954 

3501 

5281 

• 

9 

0.6096 

2731 

2371 

3731 

• 

• 

0.6096 

3456 

2938 

4549 

• 

• 

1.0668 

4849 

4039 

6310 

9 

• 

1.0668 

4012 

3229 

3130 

9 

• 

1.0668 

4931 

4293 

6353 

9 

9 

1.0668 

4791 

3730 

6084 

• 

9 

1.524 

5983 

4733 

7642 

• 

9 

1.324 

5053 

3901 

6384 

• 

T 

1.524 

5218 

4199 

6698 

• 

3.048 

1.524 

3191 

4035 

6575 

EIEIHi 

0.1324 

689 

716 

994 

• 

i 

0.1524 

783 

792 

1113 

• 

9 

0.1324 

676 

689 

966 

• 

9 

EIQIIH 

774 

732 

1079 

9 

9 

0.6096 

EQEHi 

2811 

4359 

T 

t 

0.6096 

3127 

2736 

4135 

0.1016 

2.134 

0.6096 

— 

— 

— 
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O  o  ol  olo  GO  00  BO  O  to  a  CO  O  GO  to  a  GO  op  O  OP  opppppppppppp 


Floe  Floe  Total 

thick.  length  Vel.  force 

Dam.  (m)  (m)  (m/s)  (kN) 


Damage  codes  in  damage  column: 

O  -  No  visible  damage 

C  -  Beam  cracked 

B  -  Beam  broke 

Floe 

Floe 

Floe 

Floe 

Total 

thick. 

length  Velocity  X-force  Z-force 

Total 

thick. 

length 

Vel. 

force 

(m) 

(m)  (m/s)  (N)  (N) 

(N) 

Dam.  (m) 

(m) 

(m/s) 

(kN) 
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Damage  codes  in  damage  column: 


O  -  No  visible  damage 
C  -  Beam  cracked 
B  -  Beam  broke 


Floe  Floe  Floe  Floe  Total 

thick,  length  Velocity  X-force  Z-force  Total  thick.  length  Vel.  force 

(m)  (m)  (m/s)  (N)  (N)  (N)  Dam.  (m)  (m)  (m/s)  (kN) 


Damage  codes  in  damage  column: 


O  -  No  visible  damage 
C  -  Beam  cracked 
B  -  Beam  broke 

Floe  Floe 

thick,  length  Velocity  X-force  Z-force  Total 
(m)  (m)  (m/s)  (N)  (N)  (N) 


1 0.0308 

1.219 

0.6096 

1072 

921 

1413 

mm 

i 

1.0668 

1174 

894 

1476 

• 

• 

1.0668 

1246 

872 

1320 

• 

• 

1379 

912 

1633 

• 

• 

1.0668 

1334 

988 

1660 

• 

• 

1.324 

1613 

1023 

1912 

• 

• 

1.324 

1481 

796 

1682 

• 

m 

1.324 

1668 

930 

1910 

• 

1.219 

1.324 

1384 

934 

1839 

• 

0JO3 

0.1324 

31 

36 

47 

• 

i 

0.1324 

40 

36 

34 

• 

• 

0.1324 

49 

36 

61 

• 

• 

0.1324 

40 

36 

34 

• 

• 

0.6096 

147 

163 

221 

• 

• 

0.6096 

173 

142 

224 

• 

• 

0.6096 

142 

163 

218 

• 

• 

0.6096 

178 

120 

213 

• 

• 

1.0668 

218 

222 

311 

• 

• 

1.0668 

222 

191 

293 

• 

• 

1.0668 

231 

218 

318 

• 

• 

1.0668 

231 

227 

324 

• 

• 

U24 

302 

178 

351 

• 

• 

1JI24 

400 

307 

304 

T 

T 

1.324 

378 

271 

465 

0.0308 

0J05 

1.324 

338 

347 

484 

0.0254 

3.048 

0.1324 

267 

283 

390 

ni 

m 

0.1324 

316 

298 

434 

• 

• 

EnZH 

227 

249 

337 

• 

• 

0.1324 

320 

323 

436 

• 

• 

0.6096 

347 

363 

303 

• 

• 

0.6096 

_ 

• 

• 

0.6096 

336 

343 

494 

• 

• 

0.6096 

480 

316 

705 

• 

• 

1.0668 

707 

303 

868 

• 

• 

1.0668 

347 

434 

711 

• 

• 

634 

467 

804 

• 

• 

K^Mi 

369 

414 

704 

• 

• 

1.324 

883 

494 

1014 

T 

t 

1.324 

649 

436 

782 

0.0234 

3.048 

1.324 

943 

627 

1133 
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03  00  CO 


Total 

force 

(kN) 


Damage  codes  in  damage  column: 


O  -  No  visible  damage 
C  -  Beam  cracked 
B  -  Beam  broke 

Floe  Floe 

thick.  length  Velocity  X-force  Z-force  Total 
(m)  (m)  (m/s)  (N)  (N)  (N) 


0.0254 

3.048 

1.524 

792 

427 

900 

0.0254 

QEm 

0  ’524 

258 

258 

365 

i 

0.1524 

258 

240 

353 

• 

• 

0.1524 

231 

258 

347 

• 

• 

0.1524 

254 

240 

349 

• 

• 

0.6096 

632 

694 

938 

• 

• 

0.6096 

— 

_ 

_ 

• 

• 

0.6096 

431 

418 

601 

• 

• 

0.6096 

338 

360 

494 

• 

• 

1.0668 

587 

601 

840 

• 

• 

1.0668 

383 

276 

472 

• 

• 

1.0668 

454 

320 

555 

• 

• 

1.0668 

427 

347 

550 

• 

• 

1.524 

881 

512 

1019 

• 

• 

1.524 

743 

543 

920 

• 

m 

1.524 

676 

525 

856 

• 

2.134 

1.524 

810 

725 

1087 

• 

1.219 

0.1524 

187 

200 

274 

• 

i 

0.1524 

173 

178 

249 

• 

• 

0.1524 

156 

169 

230 

• 

• 

0.1524 

182 

214 

281 

• 

• 

0.6096 

449 

534 

698 

• 

• 

0.6096 

365 

343 

500 

• 

• 

0.6096 

369 

396 

541 

• 

• 

0.6096 

338 

302 

454 

• 

• 

1.0668 

525 

472 

706 

• 

• 

1.0668 

596 

485 

768 

• 

• 

516 

387 

645 

• 

• 

1.0668 

436 

338 

552 

• 

• 

1.524 

636 

387 

745 

• 

• 

1.524 

636 

351 

727 

• 

T 

1.524 

498 

391 

634 

• 

1.219 

1.524 

667 

347 

752 

• 

0.305 

0.1524 

31 

49 

58 

• 

mm 

0.1524 

49 

62 

79 

• 

• 

0.1524 

40 

62 

74 

• 

• 

0.1524 

40 

49 

63 

• 

• 

0.60^ 

125 

160 

203 

T 

T 

0.6096 

142 

151 

208 

0.0254 

0.305 

0.6096 

125 

129 

179 

oa  n 


Floe  Floe  Total 

thick.  length  Vel.  force 

Dam.  (m)  (m)  (m/s)  (kN) 


Damage  codes  in  damage  column: 

O  -  No  visible  dcunage 
C  -  Beam  cracked 
B- Beam  broke 

Roe  Total 

length  Vel.  force 

(m)  (m/s)  (kN) 


Floe  Floe  Floe 

thick,  length  Velocity  X-force  Z-force  Total  thick, 

(m)  (m)  (m/s)  (N)  (N)  (N)  Dam.  (m) 


0.0254 

0.305 

0.6096 

125 

129 

179 

B 

0.127 

1.53 

1.36 

22.4 

i 

1 

1.0668 

182 

187 

261 

B 

i 

i 

mm 

32.6 

• 

• 

1.0668 

191 

173 

258 

B 

• 

• 

2.39 

32.3 

• 

• 

1.0668 

214 

276 

349 

B 

• 

• 

2.39 

43.6 

• 

• 

1.0668 

187 

200 

274 

B 

• 

• 

2.39 

34.2 

• 

• 

1.524 

182 

305 

B 

• 

• 

3.41 

38.1 

• 

• 

1.524 

258 

173 

311 

B 

• 

• 

3.41 

38.9 

T 

T 

1.524 

261 

205 

336 

B 

T 

_I _ 

3.41 

42.0 

0.0254 

0.305 

249 

200 

320 

B 

0.127 

1.53 

3.41 

39.9 
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APPENDIX  C:  PROCEDURE  FLOW  CHART 


Experiments 

Floating  beam  experiments 
Vertically  applied  force  at  one  end 
Beam  displacements  measured 


Numerical  Model 

Fluid  influence  coefficient  matrix 
developed  using  wave  equation  and 
acoustic  pressures  (McGilvary  et  aL) 
Ice  designated  as  elastic  material 


Numerical  predictions 

compared  with  - 

experimental  results 


Ice  impact  experiments  performed 
All  combinations  of  4  ice  lengths, 
thicknesses  and  velocities 


Comparisons  acceptable 


ABAQUS  finite  element  analysis 
used 

Ice  is  elastic  beam,  bridge  is  45° 
surface 

with  springs  for  stiffness.  Dymamic 
fluid 

foundation  determined  above. 


Numerical  predictions 

cooqjared  witti  _ 

experimental  results 


Ice  and  structure  b^avior 
assessed  using  comparisons 


Conclusions 
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